Interface depinning versus absorbing-state phase transitions 
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According to recent numerical results from lattice models, the critical exponents of systems 
with many absorbing states and an order parameter coupled to a non-diffusive conserved field 
coincide with those of the linear interface depinning model within computational accuracy. In this 
paper the connection between absorbing state phase transitions and interface pinning in quenched 
disordered media is investigated. For that, we present a mapping of the interface dynamics in a 
disordered medium into a Langevin equation for the active-site density and show that a Reggeon- 
field-theory like description, coupled to an additional non-diffusive conserved field, appears rather 
naturally. Reciprocally, we construct a mapping from a discrete model belonging in the absorbing 
state with-a-conserved-field class to a discrete interface equation, and show how a quenched disorder 
is originated. We discuss the character of the possible noise terms in both representations, and 
overview the critical exponent relations. Evidence is provided that, at least for dimensions larger 
that one, both universality classes are just two different representations of the same underlying 
physics. 
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I. INTRODUCTION 

Phase transitions separating a non-trivial from a frozen 
phase, in which the dynamics is completely arrested, ap- 
pear in a large variety of situations in physics, as well as 
in many other disciplines A central problem from 

a theoretical viewpoint is to understand how the symme- 
tries and conservation laws of the dynamics are reflected 
in the categorization of models into universality classes. 
There are two main general contexts in which this type 
of frozen states appear: 

(i) Lattice models with discrete particles; typically par- 
ticles originate "activity" and the frozen state, without 
activity is referred to as "absorbing state" This 
group appears in various disguises as cellular automata 
pj , reaction-diffusion systems , directed-percolation- 
type models || , or the fixed energy ensemble of sandpilc 
cellular automata 0], among many other examples. 

(ii) Elastic interfaces in random environments. In this 
second group, the dynamics is frozen whenever the inter- 
face is pinned by the disorder, while the non-trivial phase 
is the moving or depinned one 

The number of physical realizations of both of these 
two generic families of phase transitions is huge Jl]-||,[|(7| • 

The most prototypical universality class in the first 
group is that embracing, among many other models and 
systems, directed percolation (DP) At a contin- 

uous level the DP class is represented by the Reggeon 
Field Theory (RFT) §J, which can be written in terms 
of the following Langevin equation: 



d t p(x, t) = ap- bp 2 



w 



(Ty/p T)(x.,t) 



(1) 



the minimal field theory capturing the relevant ingredi- 
ents of the DP universality class. It can be renormalized 
using standard field theoretical methods and the associ- 
ated critical exponents can be computed in e-expansion 
||. Other universality classes of absorbing-state phase 
transitions have been identified; all of them owe their ex- 
istence to the presence of some additional symmetry or 
conservation law. Among them some example are: the 
conserved parity (CP) class, in which there are two Z^- 
symmetric equivalent absorbing states dynamical 
percolation 10 1, and the different classes of transitions 
with extra conservation laws |TT- 13|. 

In the group of pinned interfaces, the simplest con- 
tinuous model for depinning is the quenched Edwards- 
Wilkinson (QEW) equation, also called, "Linear Inter- 
face Model" (LIM) W 



(2) 



d t h(x., t) = iA7 2 /i(x, t) + F + t/(x, h) 



where p is an activity field, a, b, and a are constants and 
ij is a delta-correlated Gaussian white noise. The RFT is 



that describes an elastic interface (the Laplacian) at the 
reference height ft,(x, f), with surface tension v. under 
the influence of a constant external driving term F, and 
a quenched noise r\. Equation (^) exhibits a depinning 
transition at a critical force F c ; the interface configu- 
ration and dynamics develop critical correlations in the 
vicinity of the critical point. The standard approach for a 
theoretical analysis of the LIM is the functional renormal- 
ization group method. One-loop expressions for the mini- 
mal set of exponents have been computed by Nattermann 
et al. jlij on one hand, and by Narayan and Fisher 
(see also the more recent work by Le Doussal and collabo- 
rators Jl6|]). Here one enters technically and conceptually 
difficult terrain due to the renormalization of the whole 
disorder correlator. The outcome is that for noise fields 77 
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which do not exhibit extra translational symmetries, the 
expected depinning behavior follows, very generally, that 
resulting from a random-field uncorrelated noise term: 
the LIM universality class (l|jl|Jl|Jl||. Other univer- 
sality classes in the interfaces- in-random-media realm are 
the quenched Kardar-Parisi-Zhang equation and 
the Edwards- Wilkinson equation with columnar noise 

Recent investigations (motivated by the analysis of 
sandpile models (2^J|^] , the archetype of systems exhibit- 
ing self-organized criticality (SOC) p3|]) have demon- 
strated that different models showing a continuous transi- 
tion into an absorbing phase and with an order parameter 
coupled linearly to an extra, non-diffusive conserved field 
(NDCF) belong to a unique universality class |l3| , p^ , p5| , 
that we will refer to as NDCF class. This class differs 
from the extremely robust DP class owing to the pres- 
ence of an additional conservation law J12| . Moreover, 
the critical exponents of this class seem, within numcri 
cal accuracy, equal to those of the LIM class P, |l3| , |2"4| , p5 
This might be surprising at first sight, as in these mode 
there is no quenched disorder, as there is in LIM, and 
disorder is usually a relevant perturbation when it comes 
to universality issues. 

From a different perspective this observation is not so 
surprising, as different often tentative results have been 
reported in the literature in order to relate the dynamics 
of sandpiles to that of elastic manifolds in random media 
p6| , p7| |. Furthermore, there is one more viewpoint from 
which the coincidence between both types of models is 
not so striking, namely, that provided by the "Run-Time 
Statistics" technique |2qj . This technique or theory es- 
tablishes that quenched disorder can be mapped rather 
generically into long-range temporal correlations (i.e. a 
long-term memory) in the activity field, (note the idea 
works also the other way around) |3CJ], and has been re- 
cently applied with success to the Bak-Sneppen model 
among others p9[ . In the NDCF class the presence of 
a conserved field plays the role of a long-memory term, 
and therefore it comes not as a big surprise that it is 
equivalent to a class with quenched disorder. 

In this article we discuss in detail the relation between 
the two presented groups of transitions, i.e. absorb- 
ing states with a conserved field and pinned interfaces 
in random media, including annealed (or thermal) and 
quenched disorder respectively. The connection between 
absorbing state models in the DP class (without a con- 
served field) and their interface representation has also 
been recently considered in the literature |Q . In partic- 
ular the RFT was mapped into rather unusual interface 
equation, not resembling any known intcrfacial problem. 

The paper is structured as follows: We start in Sec- 
tion [l| by presenting the RFT-like Langevin equation for 
the recently introduced NCDF class. In III we present a 
prototypical interface model in the LIM class, in partic- 
ular the cellular automaton by Leschhorn [jljj (see also 
p8[) and work out a derivation of a Langevin equation 
for the activity density p, paying particular attention to 



the way by which the noise can be found. In Section 
IV we proceed conversely: we employ a discrete map- 
ping of a model with absorbing states in the NDCF class 
into a continuous interface representation. We end up 
with an interface equation, with several quenched noise 
terms that reflect the microscopic rules and the thermal 
noise applied in them. We discuss at this point the noise 
correlations that arise and their relevance, with the aid 
of the renormalization group (RG) literature. Finally, we 
present a discussion and an appendix in which we outline 
the relations between the exponents in the two different 
pictures. 



II. THE NCDF FIELD THEORY 

One particular system in the NDCF class (out of the 
many studied [ fL3|j2~| |) is a two-species reaction-diffusion 
model, in which one of the species in immobile (see 
. IV for a detailed definition) . It has the great ad- 
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vantage of allowing for a rigorous derivation of a coarse- 
grained field theory (or, equivalently, a Langevin equa- 
tioii) via a Fock space representation of the dynamics 
p2| , pd 24 1 . The result is in the form of a Reggeon field 
theory coupled to an extra conserved non-diffusive field, 
or what is equivalent, a RFT equation with an extra 
non-Markovian term |2^^,^|. Quite remarkably this 
Langevin equation coincides (up to irrelevant terms) with 
the one proposed previously, based only on symmetry and 
relevancy arguments, as the minimal Langevin equation 
capturing the physics of NDCF, namely f>|,p4j: 



dp(x, t) 
dt 



= ap(x) - bp(x.) 2 + V 2 p(x, t) - pi/j(x, t)p(x, t) 



£V 2 p(x,i), 



(3) 



^(x,t) 
dt 

plus higher order terms, irrelevant from naive power 
counting analysis Q. Note that the second equation, de- 
scribing the evolution of the background conserved field 
(coarse grained representation of the total number of par- 
ticles, which is conserved in the microscopic model), rep- 
resents an static non-diffusive field: in the absence of 
activity its dynamics is frozen. Observe also that the 
second equation, being linear, can be integrated out, and 
a closed equation for the activity written down. More 
concretely 



ip(x,t) = ?A(x,0) + L> f dt'^p^t 1 ). 
Jo 



(4) 



The first contribution in Eq([|), a quenched (columnar) 
disorder, represents the initial condition, while the sec- 
ond is a non-Markovian term. The Langevin equation 
(^), even though it looks rather similar to the RFT, has 
resisted all renormalization attempts; therefore predic- 
tions about critical exponents coming from an epsilon 
expansion calculation are not available so far. 
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III. PHENOMENOLOGICAL ACTIVITY 
DESCRIPTION OF LIM MODELS 

We consider a representative of the LIM class, namely 
the Leschhorn-Tang (LT) cellular automaton jlTj . In or- 
der to study its relation with standard systems with ab- 
sorbing states, we intend to cast it into a Langevin equa- 
tion describing the evolution of an activity field [|| . 

The LT automaton is defined as follows. The inter- 
face field /i(x) satisfies at each discrete time step ij the 
following equation: 



h(x, t 



i+l, 



_ (h(x,U) + l, /(x,ii)>0 



/(M,)<0 



(5) 



where the force / is given by the combination of elasticity 
and a random quenched pinning force as 



/(Mi) = V 2 /i(x,i 4 ) +??(x,/i) 



(6) 



where V 2 /i(x) is the discrete Laplacian, i. e. J2 nn h(nn) — 
2Dh(x) where nn denotes the nearest neighbors on a 
hyper-cubic lattice. A reasonable choice for the noise 
is 



r?(M) = { 



+1, 
-1. 



P 

1-p 



(7) 



when p is a random number uniformly distributed be- 
tween zero and unity. This choice implies that the aver- 
age driving force is F = (/) = 2p — 1. F plays the role 
of a control parameter. The critical point is estimated to 
be at p c ~ 0.800 @. 

At every time step, and at each site where the total 
driving force exceeds its threshold value, i.e., at each 
interface-site advance, we define an activity variable and 
set it equal to one. On the other hand, in the remaining 
lattice sites the corresponding activity takes a zero value. 
Additionally, we also define at each site and time, a con- 
tinuous "background" variable, equal to V 2 /i(x, t) + F. 
This controls the probability of each interface site to ad- 
vance at each time, regardless of whether it actually slips 
or not. Let us emphasize that this background variable 
is a conserved magnitude, i.e., it takes a constant value, 
equal to F when integrated (summed) over the whole lat- 
tice. However, locally, it favors or inhibits the generation 
of new activity. We now build up a couple of equations 
for the evolution of the two fields: the activity, p(x,t), 
and the background field, tp(x.,t), which are the coarse 
grained field analogous of the previously defined site 
variables. Using the identification between activity and 
ready-to-advance sites: h(x,t) = J Q dt'p(x,t') + h(x, 0). 
Let us write down a couple of mean-field equations for 
the two defined fields: 



dp(x, t) 
dt 

ip(x,t) = W 2 h(x,t)+F 



p(x, t) + p(x, t) G [^(x, t)] V 2 p(x, t) (8) 
2 /i(x,i) +F 

[ dtV 2 p(x,t) + V 2 /i(x,0) + F (9) 
Jo 



The justification of the different terms is as follows: 

• The term — p(x, t) describes the decay of ac- 
tive sites, that after the corresponding interface- 
advance become, in general, non-active. At a coarse 
grained level higher order corrections, as — 6p 2 (x, t) 
may also appear. In particular, they might play an 
important role in order to prevent the activity from 
growing unboundedly, i.e. in stabilizing the theory. 

• +p(x.,t)G[ijj(x,tj\ represents the fact that activity 
is created in regions where some activity is already 
present, and the rate of creation at each point is 
a function of the local background field, ip(x,t). 
Observe that the total contribution of this term 
when integrated over the whole space has to be 
zero, but locally it fosters or inhibits the creation 
of further activity. Again, higher order powers of 
yo(x, t) might also be included. 

• V 2 /?(x, t) describes the diffusion of activity. This 
terms appears generically for diffusive systems at a 
coarse grained scale. 

• In what respects the ip(x.,t) field, Eq.(^), we have 
just written its definition by equating /i(x, t) to the 
number of "topplings" (or activity events) at that 
point in all the preceding history, plus its initial 
value. 

Expanding Q[tp(pc, t)] in power series, and keeping 
only the leading contribution, we are left with a term 
+Xphc, i)ip(x, t) (where A is a constant) on the r.h.s. of 
Eq. (||) (observe that the constant term in the Taylor ex- 
pansion has to be zero as its integral has to be conserved, 
as argued before). A posteriori, we shall show that the 
omitted terms, as well as higher order corrections to 
the Laplacian term, are irrelevant in what respects large 
scale, asymptotic, properties. 

In order to account for the system fluctuations (com- 
pletely ignored so far) we now introduce a noise field 
contribution to Eq.(||). For that, as it is well known 
in field theoretical descriptions of systems with absorb- 
ing states [0|§, a RFT noise term: <r^/p(x,t)ry(x, t) is 
needed, where er is a constant and r/ a Gaussian white 
noise. This just reflects the fact that, as p is a local coarse 
grained variable its local fluctuations are proportional to 
its square-root (see [[[]-[| and references therein). It also 
captures the physical key ingredient: wherever activity 
vanishes locally, fluctuations are canceled 

Before proceeding further, let us now discuss why the 
quenched disorder of the microscopic model can be repre- 
sented by an annealed noise in our description. The key 
point is the observation that in active regions, i.e. where 
the interface advances, a new noise variable is selected 
at every time step and, as the interface does not return 
to already passed regions, there is no need to store the 
microscopic noise history, and the noise can be freshly 
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extracted from its probability distribution after every in- 
terface advance. Therefore, it is rather obvious that in 
depinned (active) regions, quenched and annealed noises 
are fully equivalent. More subtle is the connection of the 
two types of noises in what respect pinned (absorbing) 
regions. While the annealed noise, rj changes in time 
even if there is no activity in a given region, its vari- 
ations are completely irrelevant as the noise amplitude 
appears multiplied by yfp = 0. Noise (including its activ- 
ity dependent amplitude) at a given spatial point changes 
only whenever activity arrives to it, mimicking perfectly 
what happens in the microscopic interface model, where 
pinned regions can be depinned only under the presence 
of neighboring moving regions. Therefore, the consid- 
ered time-dependent noise, reproduces properly all the 
properties of the original quenched disorder. 

All previous considerations lead finally to the following 
Langevin equation for the activity field: 



dp(x, t) 
dt 



= [-1 + XF + AV 2 /i(x, 0)]p(x, t) + V 2 p(x, t) 



+Ap(x,t) / dtV 2 p{x, t') + ayjp(x, i)?7(x, t) (10) 
Jo 

where we have substituted ip by its expression com- 
ing from Eq.(^J). In general, the system is expected to 
lose memory of the initial state for long enough times, 
therefore the dependence on V 2 /i(x,0) is expected to be 
washed out. However, in some cases, as for instance one- 
dimensional systems, due to the meager phase space, and 
the slow relaxation of the initial condition, this might not 
be the case Q. 

Performing a perturbative, diagrammatic study of the 
previous Langevin equation it is easy to see (already at 
one loop level) that a new non-linearity (vertex), with the 
same degree of relevancy as the nonlinear terms already 
present in the theory (i.e. the non-local-in-time vertex 
and the noise one) is perturbatively generated: p 2 (x, i). 
In fact, this term could have been introduced also at a 
mean field level, as pointed out before, as a stabilizing 
term for the activity equation. 

Including all the discussed terms into the equation for 
p, and integrating the equation for tp, we finally obtain: 

9p{ Q^ = -<V( X > *) - t) 2 + Ap(x, t) J dt'pix, t') 

+AV 2 /i(x, 0) + V 2 p(x, t) + a v / p(x,t)r]{x, t) (11) 

where a = — 1 + FA and b > are constants. At this 
point, it is a rather straightforward exercise to verify 
that no further relevant terms are generated when in- 
cluding perturbative (diagrammatic) corrections to the 
bare theory. Therefore, the resulting Langevin equation is 
identical to the one proposed for systems with an infinite 
number of absorbing states and an activity field coupled 
to an static conserved field Eq. f3|,[i"3 24 1 . 

Summing up, we have mapped a microscopic model be- 
longing in the LIM class to the Langevin equation char- 
acterizing the NCDF class. Though our derivation is not 



rigorous, we believe it provides strong evidence that in 
fact LIM and NDCF define the same universality class. 



IV. MAPPING A REACTION-DIFFUSION 
MODEL TO DEPINNING 



In this section we proceed conversely to the previous 
one: starting from a microscopic model in the NDCF 
class we map it onto the LIM continuous equation, 
Eq.(||). To that end we follow a recipe already applied 
to many sandpile models exhibiting SOC [^tJ. Following 
p4j we consider a two-species reaction-diffusion process 
on a L d lattice, with particles of types A and B involved. 
At each site i, and at each (discrete) time step the fol- 
lowing reactions take place: 



B% > B nn , r^i — 
Ai + Bi -»• 2 x Bi , ri 
Bi -> Ai , r%. 



(12) 
(13) 
(14) 



The Ai, Bi denote particles of each kind at site i. the 
r's are the probabilities for the microscopic processes to 
occur: diffusion, rd, activation r±, and passivation r%. 
Without loss of generality we will fix = 1, implying 
that, after having the chance to react, B particles dif- 
fuse with probability one. Thus one can define a phase 
boundary between the active and absorbing phases in 
terms of the r\, r-i probabilities, with a phase transition 
in-between. We assign occupation numbers tia,%i TiB,i to 
each site. As the A particles are non-diffusive, this sys- 
tem has an infinite amount of absorbing states defined 
by nB,i = for all i, with ua,i arbitrary. 

Now we define (analogously to what is done for sand- 
piles (2?]||) a height field H(x,t) which increases by one 
unit every time a site gives one (or more than one) ac- 
tive, diffusing B particle to one (or more than one) of its 
neighbors. When this happens we say, using the sand- 
pile terminology, that the site "topples" . In this way, the 
ii-field measures the integrated activity at x up to time 
t. 

The mapping to an interface automaton with quenched 
noise is based on the fact that both, the reactions between 
the A and B species, and the diffusion of particles can 
be accounted for by looking at their net effects at every 
time B particles leave the site x. One just have to look at 
ha and ns when the site becomes active and a particle 
diffuses out. The dynamics of H can be written as: 



1 ' + j l#(x,i), /(x,F)<0 



(15) 



which is formally identical to the Leschhorn automaton 
in the LIM class, with a local "force" defined as 



/(x,J7)=n tot (x,JT)-£(x,J?) 



(16) 
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where n to t(x, H) = n^(x, i?)+ns(x, H) is the total num- 
ber of particles at x, and £ is a local random threshold 
that results from the microscopic processes. More con- 
cretely, the noise £ is defined as follows: Consider the 
site x after the H-th toppling, either ^(x, H) = or 
n b (x, H ) > (this last can be the case if and only if 
particles have arrived from the nearest neighbors at the 
same time step). In the first case, it will remain zero until 
a particle arrives from a nearest neighbor site; then one is 
free to choose a value for £ (x, H ) such that it makes the 
force / negative in the time interval between topplings H 
and H+l. In the second case, ns(x) will fluctuate owing 
to the microscopic passivation and activation processes, 
either going to iis(x) = or inducing a toppling at the 
next time step. The relative probabilities of these two 
alternatives, as derived from the microscopic dynamics, 
are captured in the £(x, H) probability distribution. 

Observe that £ depends solely upon the total number of 
particles after the preceding toppling and the microscopic 
dynamical rules. In particular, the larger n to tai the larger 
the probability to have many B particles and the larger 
the probability to topple. Let us also remark that the im- 
mobile grains ua constitute a "pinning force" (the larger 
their relative number, the lesser the probability to top- 
ple) . The point- wise noise field £ (x, H) should have weak 
two-point correlations in x since, in particular, it is de- 
pendent on the number of grains received from the nn's at 
the interface location H(x) which induces weak site-site 
correlations. The fact that ua changes slowly will make 
the ff-part of the noise correlator (£(x, H)£(x , H )) less 
trivial than a simple delta- function S(H — H'). 

Equation ([l^) can be considered as a discrete interface 
equation 

A H 

-^ = 9(f(x,t)). (17) 

It can be rewritten with the help of two particle-fluxes: 
n™ and n™*, are the number of grains added to or re- 
moved from a given site x up to time t, respectively. 
Let us also define g as the average number of particles 
given to the nearest neighbors at each toppling event. It 
is clear that for long enough times n° Mt w gH(yi); rela- 
tive deviations from this equality being negligible asymp- 
totically. Defining also the average value of n™, n™ , 
as = <7/2dX)x„„ H( x nn,t), we can compute a noise 
t(x, H) as the deviation of with respect to its average 
value: 

r(x,*)=<-^£ff(x nn ,*). (18) 

In other words, r counts the relative proportion of parti- 
cles diffused out from the neighbors that actually arrive 
to the site under consideration, compared with its av- 
erage value. A site to which particles have toppled in 
excess will take a positive value of r, and therefore will 
be more likely to topple in the following time steps. 

Plugging this into Eq. (|l6|) , and using that rHot(x, t) — 
n tot (x, 0) + ri™ — ri° M * , we can write [gj 



f =2d V 2 H + F( X ,0)-ax,H)+T(x,H) (19) 

where F(x, 0) = n tot (x, t = 0). 

The discretization in Eq.(|l7|) can be understood so 
that the rules result in an effective force /' that is ex- 
actly unity when the interface field H advances. Thus 
AH/ At ee f'9(f) = f'e(f) |3- This construction can 
be achieved by picking £ to have exactly the right value 
in order to make the force driving the interface equal to 
unity, if it is larger than zero. One arrives finally at the 
discretized interface equation 

= 2d + F(X ' 0) " $(X ' H) + T(X ' H) - (20) 
Let us stress the presence of three different noise terms: 

1. F(x,0) represents the original total- number-of- 
particle configuration at t = 0, and is therefore 
a columnar noise term [ ^o| . It induces an ini- 
tial transient regime until eventually, the dynam- 
ics washes out the dependence of the original con- 
figuration. In general, columnar disorder is irrel- 
evant in the renormalization group sense as com- 
pared to quenched noise; therefore using relevancy 
arguments, it could be eliminated, at least in high 
enough dimensions, close or above the critical one 
d c = 4. Notice that this statement is equivalent 
to the LIM symmetry, by which static force fields 
F(x, 0) (independent of H) is completely equivalent 
to the existence of a non-trivial initial interface pro- 
file H(x,t = 0). However, in low dimensional sys- 
tems, and in particular in d = 1, due to the meager 
phase space, relaxation times might be huge, and 
the time needed to eliminate the dependence on the 
initial particle distribution divergently large [ p4[ . 

2. The noise term £(x, H) represents the local thresh- 
old, determining whether a site with some B parti- 
cles topples at a given time or, alternatively, they 
are transformed into A particles by microscopic 
processes. It captures the in-site microscopic dy- 
namics, and depends essentially on n tot , and on the 
microscopic probabilities. On a nutshell, it says 
how many of the n to t particles are of type A af- 
ter the microscopic dynamics has operated in the 
corresponding time step: if all n to t are of type A 
then £ > n to t, and / < 0; conversely, if any of the 
particles is of type B then £ < n to t and / > 0. Ob- 
serve that if the diffusion probability was smaller 
than unity, then we should substitute £(x, H) by 
a "thermal noise" £(x, t), i.e. £ would change its 
value after every time step instead of changing only 
after each toppling: this is due to the fact that if 
rd < 1 then a site x including B-type particles could 
not topple at time t (£(x, t) below threshold), and 
do so at a future time t' (£(x, t') above threshold). 
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This "thermal noise" would generate a rounding off 
of the transition, but the critical exponents should 
not be affected by this irrelevant perturbation Q . 
Therefore, we stick to the simplest case, r& = 1. 

3. The noise term r keeps track of the Brownian mo- 
tion of particles; i.e. it takes into account the fact 
that particles are not homogeneously distributed 
among the nn, but one of them is picked up ran- 
domly for each toppling event. It changes slowly 
since the effect of the random choices (directions) 
on the configuration is slow. This is in particular 
true since the noise r is conserving, as the num- 
ber of particles is conserved (and as can be seen by 
integrating Eq. (|l8|)). A key point is that, analo- 
gously to what discussed in the preceeding section, 
the choice to give a particle to a certain neighbor 
can be taken to be "quenched" , i. e. chosen in ad- 
vance at t = 0, or "annealed" , i. e. decided on the 
spot. The correlator of r can be generically written 
as 

(r(x , ff'Mx, H)) ~ /|| (x - x)/ x (#' - H). (21) 

The (so far unknown) correlators f\\ and f± reflect 
the discrete nature of the choices in the dynam- 
ics. Two microscopic reasons lead immediately to 
non-trivial correlations in r: 

(a) The noises r at the nn's of site x are correlated 
due to an exclusion effect: If a site gives out a dif- 
fusing B particle to a neighbor, then all the other 
neighbors are excluded. The actual coarse-grained 
noise correlations are harder to assess, since the 
fluctuations in the particle flux that r measures 
make the interface to fluctuate, and thus a sepa- 
rable noise correlator as Eq. (|2l]) is hard to com- 
pute. The easiest way to analyze the correlations 
among the different sites is therefore to compute 
the noise correlator from numerics of the micro- 



scopic model, using the noise definition Eq. (18) 



This programme has been pursued for sandpiles 
@- 

(b) At each site the noise follows the dynamics of a 
random walker. In fact, every time a nearest neigh- 
bor topples, the choice (give the particle to x or to 
a different site) makes it so that f± ~ (H — H) 1 / 2 
since at every step r can go "up" or "down" with 
respect to the average. 

Therefore, reciprocally to what done in the previous 
section, we have mapped the reaction-diffusion process 
into an interface equation. The dynamics of this inter- 
face equation follows exactly the history of a reaction dif- 
fusion process, the details of which are mapped into the 
quenched noises £ and t, and a columnar noise -F(x, 0). 
Let us remark that the existence of a conservation law 
has played a key role in order to obtain a Laplacian in 
Eq.©. 



Finally, using standard renormalization group argu- 
ments about the relevancy of different operators, we 
can eliminate higher order irrelevant terms and noise- 
correlations, and then we are left with the LIM equation 
for point- disorder ]l4| , |l6| ] (see also the Appendix). 

It must be emphasized, that the mapping works in both 
ways, it is evident that the noise construction can be in- 
verted to yield a reaction diffusion process, that corre- 
sponds to an interface model, assuming that the original 
noise terms have the right correlation and conservation 
properties. The interface model Eq.(20) resembles very 
much the one that corresponds to the Manna sandpile 
automaton, with the addition of the £-noise term which 
is more point-disorder like than the r-term. 

Summing up, reciprocally to what was done in the pre- 
vious section, in this one, we have constructed a mapping 
between a microscopic model in th NCDF class into the 
Langevin equation for the LIM class. 



V. DISCUSSION 



We have presented strong theoretical evidence that, 
rather generically, the universality class of systems with 
many absorbing states and order parameter coupled to a 
non-diffusive conserved field, the NDCF class, and that 
of the linear interface model with point-disorder coincide. 
This fact, already pointed out from numerical simulations 
]l3| , p4| p| is true at least nearby the critical dimension 
d c = 4, where relevancy arguments are reliable. In low 
dimensional systems this equivalence could break down 
owing to the existence, for example, of slow decaying ini- 
tial conditions J34| . For the frozen configurations in the 
point-disorder LIM it is known that the correlations of 
the forces n(x, H) acting on the interface vanish. In the 
case of NDCF models, like the Manna sandpile, such cor- 
relations (now computed from the particle configuration 
in frozen configurations) may become non-zero: this is 
a future avenue for numerical studies, but hopefully this 
would be a irrelevant feature. 

Likewise, if one considers a noise field for the LIM 
(Eq. (0)) with non-trivial (power-law) bare correlations 
in x or H, it is unclear at this point how these should 
be reflected in the construction of a Langevin equation 
for the corresponding activity field, like Eq. ( TQ) . Corre- 
lations in the local forces (or "activity thresholds") will 
affect the way the coarse-graining works. For instance, 
due to the noise structure the pinned and still-active re- 
gions will be correlated. 

In order to establish the connection between the two 
classes we have mapped a discrete interface model into 
the Langevin equation characteristic of the NDCF class, 
and conversely mapped a discrete model in the NDCF 
into the well known Langevin equation describing the 
LIM class. In order to have a more rigorous proof, one 
should be able to map one Langevin equation into the 



G 



other, but this, being the Langevin equations coarse- 
grained representations of the microscopic models, is not 
an easy task to fulfill, and remains an open challenge. 

Let us remark that a similar problem remains open; 
namely, the rigorous connection between the Quenched 
KPZ depinning transition and directed percolation 

depinning [Q,[56| in two dimensional systems, (and to di- 
rected surfaces in higher dimensions ||37|| ) . It is clear from 
numerics, that indeed these two universality classes coin- 
cide, but a satisfactory proof is, to the best of our knowl- 
edge, still lacking. 

It was the hope, that the possibility of renormalizing 
the NDCF Langevin equation using standard RG tech- 
niques, of problem from the RFT-like equation approach, 
could shed some light on the (in principle, technically 
more difficult and obscure) functional renormalization 
group analysis of the interface equation with quenched 
noise. However, the difficulties encountered in renormal- 
izing, using standard perturbative schemes, the Langevin 
equation for NDCF ]24|,^5| are considerable; and have 
made all the attempts to renormalize the theory to fall 
through. Renormalizing the NDCF Langevin equation 
and relating the derived critical exponents to those ob- 
tained using functional RG for LIM remains an open and 
very challenging problem. 

Finally let us also point out that all the discussions pre- 
sented in this work deal with the "constant force" (in the 
interface language) or "fixed energy" (in the absorbing- 
state terminology) ensemble. They can be easily ex- 
tended to the "constant force" or "slow driving" ensemble 
P.p7| , in which the system self-organizes into its critical 
state. This point is, however, not essential since all evi- 
dence points to the fact that if two models belong to the 
the universality class, they continue to share the same 
set of critical exponents upon changing ensemble. 



APPENDIX 

The scaling of the phase transition in the absorbing- 
state representation is characterized by the exponents u± , 
z/|| , z and j3. These describe the correlations in the activ- 
ity p in the spatial and time directions, the development 
of the correlations in time, and the behavior of p above 
the critical point, respectively. One has the scaling rela- 
tion 

p(A, L) = L-P/^KiL^A) , (22) 

where A is the distance to the critical point, and 1Z is a 
scaling function with TZ(x) ~ x 13 for large x. For L ^> 
£ ~ A _I/i we expect ~p~Z ~ A' 3 (here £ is the correlation 
length). When A = we have that p^(0, L) ~ L~@/ v± . 
For A > 0, by contrast, /5a approaches a stationary value, 
while for A < it falls off as L~ d . These can be used to 
establish the numerical values of the exponents. 

In the interface representation the relevant exponents 
are v, z as above, with the convention that v = v\\ . Usu- 



ally it is assumed that the dynamics is self-affine, which 
implies that v± = x v \\ lHJzl- This defines the roughness 
exponent x that characterizes the spatial correlations of 
the interface. If "simple scaling" [gSM holds, then one 
has a unique roughness exponent and we can write for 
the interface width w 

W%L)~{%: : |<£ ( 23 ) 

using also the early-time exponent (3 W . If simple scal- 
ing holds, we have the exponent relation f3 w z = a [p5| . 
If only one timescale is present, the growth exponent 
is related to the activity time-decay exponent, 8, via 

For point-like disorder the first-loop functional renor- 
malization group result reads \ — (4 — d)/3, and z = 
2 — (4 — <f)/9 [[lij ; see the extension to second order in 
Jr|. From these, using the exponent relations, the other 
exponents follow. For rather generic bare disorder cor- 
relators the implication is that the full correlator flows 
in the renormalization to this "random field" (or point- 
disorder) fixed point function, and thus the exponents are 
the same. However, numerics in particular in ID implies 
that the real exponents are different from the one-loop 
results. This has recently been explained in terms of 
two-loop corrections, but the traditional interpretation 
has been in terms of "anomalous scaling" Jl7],[39| , mean- 
ing that as t — ► oo, the typical height difference between 
neighboring sites increases without limit. 



ACKNOWLEDGMENTS 

We acknowledge partial support from the Euro- 
pean Network contract ERBFMRXCT980183, from the 
Academy of Finland's Center of Excellence Program, and 
DGESIC (Spain) project PB97-0842. We thank R. Dick- 
man, A. Vespignani, R. Pastor-Satorras, K. B. Lauritsen, 
and S. Zapperi for useful discussions and long standing 
pleasant collaborations. 



[1] J. Marro and R. Dickman Nonequilibrium Phase Tran- 
sitions in Lattice Models (Cambridge University Press, 
Cambridge, 1999). 

[2] See, G. Grinstein and M. A. Munoz, The Statistical Me- 
chanics of Systems with Absorbing States , in "Fourth 
Granada Lectures in Computational Physics", edited 
by P. Garrido and J. Marro, Lecture Notes in Physics, 
Vol. 493 (Springer, Berlin 1997), p. 223, and references 
therein. 

[3] H. Hinrichsen, Adv. Phys. 49, 1, (2000). 
[4] E. Domani and W. Kinzel, Phys. Rev. Lett. 53, 311 
(1984). 



7 



[27] 



[28] 
[29] 

[30] 



[5] A. Vespignani, R. Dickman, M. A. Munoz, and S. Zap- 
peri, Phys. Rev. Lett. 81, 5676 (1998); Phys. Rev. E 62, [25] 
4564 (2000); R. Dickman, M. A. Munoz, A. Vespignani, 
and S. Zapperi, Braz. J. Phys. 30, 27 (2000). 

[6] T. Halpin-Healy and Y.-C. Zhang, Phys. Rep. 254, 215 

(1995) . 

[7] A. -L. Barabasi and H. E. Stanley, Fractal Concepts 
in Surface Growth, (Cambridge University Press, Cam- 
bridge, 1995). [26] 

[8] H. K. Janssen, Z. Phys. 42, 141 (1981); ibid. 58, 311 
(1985); P. Grassberger, Z. Phys. B 47, 465 (1982); Cardy 
J.L., Sugar R.L. (1980), J. Phys. A 13, L423. 

[9] P. Grassberger, F. Krause, and T. von der Twer, J. 
Phys. A 17, L105 (1984); P. Grassberger, ibid. 22, L1103 
(1989); H. Takayasu and A. Yu. Tretyakov, Phys. Rev. 
Lett. 68, 3060 (1992); I. Jensen, Phys. Rev. E 50, 3623 
(1994); N. Menyhard and G. Odor, J. Phys. A 29, 7739 

(1996) ; J. Cardy and U. C. Tauber, Phys. Rev. Lett. 77, 
4780 (1996). 

[10] J.L. Cardy, J. Phys. A 16, L709 (1983); J.L. Cardy and 

P Grassberger, J. Phys. A 18, L267 (1985). 
[11] F. van Wijland, K. Oerding, and H. J. Hilhorst, Physica 

A 251, 179 (1998). See also R. Kree, B. Schaub, and B. 

Schmittmann, Phys. Rev. A 39, 2214 (1989). 
[12] If the conserved field is diffusive, there is not an infinite 

number of absorbing states, and the phenomenology is 

completely different. See also: J. E. de Freitas, L. S. Lu- 

cena, L. R. da Silva, and H. J. Hilhorst, Phys. Rev. E 61, 

6330 (2000); J.-P. Leroy, H. J. Hilhorst, K. Oerding, and 

F. van Wijland, J. Stat. Phys. 99, 1365 (2000), and 0. [31] 
[13] M. Rossi, R. Pastor-Satorras, and A. Vespignani, Phys. 

Rev. Lett. 85, 1803 (2000). [32] 
[14] T. Nattermann, S. Stepanow, L.-H. Tang, and H. 

Leschhorn, J. Phys. (France) II 2, 1483 (1992); H. 

Leschhorn, T. Nattermann, S. Stepanow, and L.-H. Tang, [33] 

Ann. Physik 7, 1 (1997). 
[15] O. Narayan and D. S. Fisher, Phys. Rev. B 48, 7030 

(1993). 

[16] P. Chauve, P. Le Doussal, and K. J. Wiese, Phys. Rev. [34] 

Lett. 86, 1785 (2001). 
[17] H. Leschhorn, Physica A 195, 324 (1993). 
[18] H. Leschhorn and L.-H. Tang, Phys. Rev. Lett. 70, 2973 

(1993). 

[19] L.-H. Tang, M. Kardar, and D. Dhar, Phys. Rev. Lett. [35] 
74, 920 (1995). [36] 

[20] G. Parisi and L. Pietronero, Europhys. Lett. 16, 321 
(1991); Physica A 179, 16 (1991). 

[21] P. Bak, C. Tang, and K. Wiesenfeld, Phys. Rev. Lett. 

59, 381 (1987); Phys. Rev. A 38, 364 (1988). [37] 

[22] S. S. Manna, J. Phys. A 24, L363 (1992); J. Stat. Phys. 

59, 509 (1990). [38] 

[23] H. J. Jensen, Self organized criticality, (Cambridge Univ. 

Press, Cambridge, 1998). G. Grinstein, in Scale Invari- [39] 
ance, Interfaces and Nonequilibrium Dynamics, NATO 
Advanced Study Institute, Series B: Physics, vol. 344, A. 
McKane et al., Eds. (Plenum, New York, 1995). 

[24] R. Pastor-Satorras, and A. Vespignani, Phys. Rev. E 62, 



5875 (2000); See also |cond-mat/010135£ . 
M. A. Munoz, R. Dickman, R. Pastor-Satorras, A. 
Vespignani, and S. Zapperi, "Sandpiles and absorbing 
state phase transitions: recent results and open prob- 
lems," to appear in Proceedings of the 6th Granada sem- 
inar on computational physics, Ed. J. Marro and P. L. 
Garrido. (American Institute of Physics, 2001). e-print: 
cond-mat/001 14471 



O. Narayan and A. A. Middleton, Phys. Rev. B 49, 244 
(1994); M. Paczuski, S. Maslov, and P. Bak, Phys. Rev. 

E. 53, 414 (1996). M. Paczuski and S. Boettcher, Phys. 
Rev. Lett. 77, 111 (1996); L. A. N. Amaral and K. B. 
Lauritsen, Phys. Rev. E 54, R4512 (1996); ibid. 56, 231 

(1997) ; C. Tang and P. Bak, Phys. Rev. Lett. 60, 2347 
(1988); D. Cule and T. Hwa, Phys. Rev. B 57, 8235 

(1998) . 

K. B. Lauritsen and M. J. Alava, preprint cond-mat/- 
9903349; M. Alava and K. B. Lauritsen, Europhys. Lett. 
53, 569 (2001); M. J. Alava and K. B. Lauritsen, in 
preparation. 

M. Marsili, J. Stat. Phys. 77, 733 (1994). 

M. Felici, G. Caldarelli, A. Gabrielli, and L. Pietronero, 

Phys. Rev. Lett. 86, 1896 (2001). 

For example, in invasion percolation the bond threshold 
is thus related to the time the bond remains active on the 
boundary, and the quenched variable is replaced with an 
annealed one, with a memory. In sandpiles, the local dy- 
namics can be mapped to quenched disorder via almost 
exactly the inverse route. 

R. Dickman and M. A. Munoz, Phys. Rev. E 62, 7632 
(2000). 

M. Doi, J. Phys. A 9, 1465 (1976); L. Peliti, J. Phys. I 
46, 1469 (1985); B. P. Lee and J. L. Cardy J. Stat. Phys. 
80, 971 (1995). 

D. J. Amit, Field Theory, the Renormalization Group 
and Critical Phenomena, (World Scientific, Singapore, 
1992). J. Zinn- Justin, Quantum field theory and critical 
phenomena, Clarendon Press, (Oxford), 1994. 
We believe that effects like the ones described here are at 
the basis of the discrepancy between the one-dimensional 
Manna model and the 1-d LIM class. R. Dickman, M. 
Alava, M. A. Munoz, J. Peltola, A. Vespig nani, and S. 
Zapperi, p reprint 2001 ( |3ond-mat/010138l| ). 

F. Family and T. Vicsek, J. Phys. A 18, L75 (1985). 
L.-H- Tang and H. Leschhorn, Phys. Rev. A 45, R8309 
(1992); S. Buldyrev et al. Phys. Rev. A 45, R8313 (1992); 
Z. Olami, I. Procaccia, and R. Zeitak, Phys. Rev. E 49, 
1232 (1994). 

A. L. Barabasi, G. Grinstein, and M. A. Munoz, Phys. 
Rev. Lett. 76, 1481, (1996). 

R. Dickman and M. A. Munoz, Phys. Rev. E 62, 7632 
(2000). 

J. M. Lopez, Phys. Rev. Lett. 83, 4594 (1999); J. M. 
Lopez and M. A. Rodriguez, Phys. Rev. E 54, R2189 
(1996); J. M. Lopez, M. A. Rodriguez and R. Cuerno, 
Phys. Rev. E 56, 3993 (1997). 



8 



